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ABSTRACT 

Details of processes through which galaxies convert their gas into stars need to be 
studied in order to obtain a complete picture of galaxy formation. One way to tackle 
these phenomena is to relate the H I gas and the stars in galaxies. Here, we present 
dynamical properties of Damped and sub-Damped Lyman-a Systems identified in H-a 
emission with VLT/SINFONI at near infra-red wavelengths. While the DLA towards 
Q0302— 223 is found to be dispersion-dominated, the sub-DLA towards Q1009— 0026 
shows clear signatures of rotation. We use a proxy to circular velocity to estimate the 
mass of the halo in which the sub-DLA resides and find M^ a ; o =10 12 6 M Q . We also 
derive dynamical masses of these objects, and find Md !/ „=10 10 3 M Q and 10 10 9 M©. 
For one of the two systems (towards Q0302— 223), we are able to derive a stellar mass 
of M„,=10 95 M© from Spectral Energy Distribution fit. The gas fraction in this object 
is l/3 rd , comparable to similar objects at these redshifts. Our work illustrates that 
detailed studies of quasar absorbers can offer entirely new insights into our knowledge 
of the interaction between stars and the interstellar gas in galaxies. 

Key words: Galaxies: - galaxies: kinematics and dynamics - quasars: absorption 
lines - quasars: individual: Q0302-223, Q1009-0026 



1 INTRODUCTION 

Tremendous progress has been made over the last decade in 
establishing a broad cosmological framework in which galax- 
ies and large-scale structure develop hierarchically over time, 
as a result of gravitational instabilities in the density field. 
The next challenge is to understand the physical processes 
of the formation of galaxies and structures and their inter- 
actions with the medium surrounding them. Of particular 

* Based on observations collected during programmes ESO 80. A- 
0330 and 80.A-0742 at the European Southern Observatory with 
SINFONI on the 8.2 m YEPUN telescope operated at the Paranal 
Observatory, Chile, 
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importance are the processes through which these galax- 
ies accrete gas and subsequently form stars (Putman et al. 
2009). The accretion of baryonic gas is complex. Recently, 
several teams (Birnboim &: Dekel 2003, Keres et al. 2005) 
have realized that, in halos with mass < p) 11 - 5 " 12 M©, bary- 
onic accretion may not involve the traditional shock heating 
process of White & Rees (1978). Similarly, details of pro- 
cesses through which galaxies convert their gas into stars 
are still poorly understood. But the observational evidences 
for accretion are scarce. A related signature is that the to- 
tal amount of neutral gas in the Universe, Qhi, is almost 
constant over most of the cosmic time (Prochaska, Herbert- 
Fort & Wolfe 2005, Peroux et al. 2005, Noterdaeme et al. 
2009), unlike the history of the star formation rate which 
peaks around z=l (Hopkins & Beacom 2006 and references 
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therein). This shows the importance of ongoing global gas 
accretion and the conversion of atomic gas to molecular gas 
in the star formation process (Hopkins, Rao & Turnshek 
2005, Bauermeister et al. 2010). 

One way to tackle these problems is to relate the H I 
gas and the stars in galaxies. While radio observations now 
provide detailed constraints on the H I content of large sam- 
ple of galaxies (Zwaan et al. 2005), they are still limited to 
redshift z~0. Conversely, the study of quasar absorbers, the 
galaxies probed by the absorption they produce in a back- 
ground quasar spectrum, is insensitive to the redshift of the 
object (Wolfe et al. 1995). Indeed, the H I content of the 
strongest of these quasar absorbers, the so-called Damped 
Lyman-a systems (DLAs), have been measured in samples of 
several hundreds of objects from the Sloan Digital Sky Sur- 
vey (SDSS). However, studying the stellar content of these 
systems has proven very challenging until now. 

In paper I of these series (Peroux et al. 2010), we re- 
ported the detection at infra-red wavelengths of H-a emis- 
sion of a Damped Lyman-a (DLA) system at z a b s =1.009 
and one sub-Damped Lyman-a system at z a b s =0.887 with 
N(H i)>10 19 (Peroux et al. 2003). Here, we take advantage 
of 3D spectroscopy made possible by the SINFONI instru- 
ment on VLT to study the spatially resolved kinematics and 
therefore the dynamical state of these galaxies. In all the re- 
maining, we assumed a cosmology with Ho=71 km/s/Mpc, 
n M =0.27 and fi A =0.73. 



2 SINFONI DETECTIONS OF H-a EMISSION 

The observations presented here were carried out at the Eu- 
ropean Southern Observatory with the near-infrared field 
spectrometer SINFONI on Unit 4 of the Very Large Tele- 
scope. First results with that instrument reported a high 
success rate in detecting metal absorption systems traced 
by Mg II absorbers (Bouche et al 2007a) at ~1, while lower 
detection rates are found at higher-redshifts z~2 (Bouche et 
al., private communication). In Paper I (Peroux et al. 2010), 
we have reported two secure detections of the redshifted H- 
a emission of high H I column density quasar absorbers: 
a DLA with log N(H l)=20.36±0.11 at z abs = 1.009 towards 
Q0302-223 and a sub-DLA with log N(H i)=19.48±g;g| at 
Zabs=0.887 towards Q1009-0026 (see Table 1). We detect 
galaxies associated with the quasar absorbers at impact pa- 
rameters of 25 and 39 kpc away from the quasar sightlines, 
respectively. 

For field Q0302— 223 where the quasar is bright enough, 
we have used the quasar itself as a natural guide star for 
adaptive optics in order to improve the spatial resolution 
(see Table 1). The two data sets have a resulting Point 
Spread Function (PSF) of 0.7" and 1.0". Both the objects 
in our study have well-determined absorption line properties 
determined from spectra of the background quasars. Thus, 
the metallicities are well-determined for both of our target 
absorbers. An estimate of the emission metallicities is made 
based on the N n/ H-a parameter (Pettini & Pagel 2004): 
the galaxy with the smaller absorption line metallicity also 
has the smaller emission-line metallicity. The absorber to- 
ward Q0302— 223 is the more gas-rich absorber, while that 
toward Q1009— 0026 is the more metal-rich absorber. Using 
the H-a luminosity we then derived the star formation rate 



assuming the Kennicutt (1998) flux conversion corrected to 
a Chabrier (2003) initial mass function. We find low star for- 
mation rates (not corrected for dust extinction), of 1.8±0.6 
and 2.9±1.0 M /yr (see Table 1). These values of star for- 
mation rates are among the lowest that have ever been possi- 
ble to detect in quasar absorber searches with ground-based 
observations at z~l and 2. 



3 KINEMATICS 

In addition to the identification and redshift confirmation 
of the galaxy responsible for the quasar absorbers, the SIN- 
FONI data allow for a study of the dynamical properties of 
the galaxies. We extracted maps of the velocity-integrated 
line fluxes, relative velocities, and velocity dispersion from 
the reduced data cubes using the following technique. For 
each pixel, we fitted the line profile to convolved gaussian 
with a template of the instrument profile following e.g., 
Forster-Schreiber et al. (2006, 2009). Thus, the fitting takes 
into account the intrinsic emission line profile. In addition, 
weighted fits are performed using the MC resampling based 
on the spectral noise profile. A boxcar median smooth of 
size 2x2 pixels was applied prior to the line fitting. The ve- 
locity resolution is 150 km/s while the resulting accuracy of 
the centroid of the line can be derived to 0.1 of a pixel (29 
km/s). These maps are displayed in Figure 1. 

Thanks to these, critical information on the dynami- 
cal state and properties of galaxies associated with quasar 
absorbers can be derived. The galaxy associated with the 
z ab s = l-009 DLA towards Q0302-223 shows little signs of ro- 
tation and significant amounts of dispersions. It has an axis 
ratio corrected for Point Spread Function (PSF) of b/a=0.47 
corresponding to an inclination of sin i=0.88. The maximum 
velocity extracted from the data cube is 2 V m ax sin i= 20 
km/s, corresponding to an inclination corrected maximum 
velocity of V ma x~ll km/s. Given that the velocity disper- 
sion (corrected from instrumental effects) is 59 km/s, the 
ratio v/<7 in this object is 0.19 (see Table 2), much smaller 
than seen in local disk galaxies (v/a— 10-20), and at high- 
redshifts (Epinat et al. 2009, Cresci et al. 2009, Lemoine- 
Busserolle et al. 2010). This might indicate that this object 
is pressure-supported, and early-type morphologies come to 
mind. 

However, the morphology of this object indicates that 
this galaxy might be a late type object. The SINFONI data 
shows that it is highly elongated. This is supported by the 
archival HST/WFPC2 images covering this field (Le Brun 
et al. 1997). In these data, two separate components can be 
seen but their redshifts cannot be estimated (see paper I 
and the two stars in Figure 1 which indicate the positions 
of the components). The two individual components are not 
resolved in the SINFONI observations and might explain the 
elongated shape or only one is seen but the detection of H- 
a emission at the position of the absorber provides secure 
identification. On interpretation is that we might be seeing 
a dispersion due to a small difference in redshift of two in- 
teracting galaxies. However, the H-a light profile appears to 
be exponential with a scale length of 0.6". Thus, alterna- 
tive explanation is that this galaxy is dispersion dominated, 
albeit it already has a disky morphology. This is indicative 
of a young flattened object, with effective radius (contain- 
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Table 1. Summary of the properties of two absorber galaxies detected with SINFONI. 



Quasar 


•Sabs 


logN(H i) 
[atoms/cm 2 ] 


[Zn/H] 


impact parameter 
[kpc] 


Av a 
[km/s] 


F(H-a) 
[crg/s/cm 2 ] 


Lum(H-a) 
[erg/s] 


SFR b 

[M /yr] 


Q0302-223 
Q1009-0026 


1.009 
0.887 


20.36+S-" 
19 48 + ' 05 


-0.51±0.12 
+0.25±0.06 
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4.1±1.4xl0 41 
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1.8±0.6 
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a : Observed velocity shift between the z^ a and z eul of the detected galaxies. 
b : The SFR estimates are not corrected for dust extinction. 
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Figure 1. H-a flux map, H-a velocity field and H-a velocity dispersion maps of the two quasar absorbers. The colour-scale 
indicates the flux in erg/s/cm 2 in the H-a flux map on the left and the velocities or velocity dispersions in km/s for the middle and 
right set of panels. North is up and east is to the left. The thin lined, black contours indicate the position of the quasar and the arrow 
represents 1 arcsec which corresponds to 8.1 kpc in the case of Q0302— 223 and 7.8 kpc in the case of Q1009— 0026. The system at the 
top (Q0302— 223) is dispersion-dominated with H-a luminosity consistent with that of a disk. The system at the bottom (Q1009— 0026) 
shows a blueshifted and redshifted component of the gas on its velocity map. This pattern is a clear characteristic of a rotating disk. 



ing half of the light) of R e =0.5" or 4 kpc, whose kinematics 
show that sustained rotation has not yet occured. In other 
words, this object will have to contract further. 

By contrast, the absorber toward Q1009— 0026 has a 
morphology and kinematics consistent with that of a disk, 
with a normal dispersion profile (v/<r=1.45) peaking at the 



center a pea k — 190 km/s and flattening out in the outer-parts 
o"disfc=60-70 km/s (see Table 2). This object is more face-on 
with an estimated inclination of sin i=0.60 derived from an 
axes ratio b/a=0.80. It shows clear signatures of rotation 
with systematic velocity gradients. Its v/a is not typical of 
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local disk galaxies which have v/<r=10-20, but the system- 
atic gradient still favors a spiral galaxy. 

We also use the kinematic maps to estimate the size of 
the systems, correcting for the seeing in quadrature. These, 
in turn, are used to compute the area of the objects assuming 
that the inclined disks appear as ellipses. Using these esti- 
mates of the sizes of the detected galaxies, we compute their 
star formation rate surface densities (see Table 2) and find 
Ss_Ffl=0.13 and 0.31 solar masses per year per kpc squared 
for absorbers towards Q0302-223 and Q1009-0026 respec- 
tively. It is interesting to note that these values are above the 
threshold at 0.1 solar masses per kpc squared where galac- 
tic winds are observed (Heckman 2002). In the log E gas -log 
SsFi? plot of Bouche et al. (2007b), our galaxies appear 
to be similar to the z~1.5 colour-selected galaxies from the 
Gemini Deep Deep Survey (Abraham et al. 2004), i.e. at 
the lower end of the relation with respect to higher-redshift 
galaxies and submillimeter-selected systems. The results of 
these are shown in Table 2. 



4 MASS ESTIMATES 

Dynamical Masses: These observations provide direct obser- 
vational estimates of dynamical masses for galaxies selected 
on neutral gas HI content. For Q0302— 223, we use the virial 
theorem to calculate the total mass of the system since the 
rotation is not established yet (Wright et al. 2009): 



M, 



dyn 



Co 2 n /a / G 



(1) 



where C is a constant (C=5 for a spherically uniform 
density profile and ri/2 is the half-light radius. We derive xi 
from a curve of growth with ellipsoidal apertures. We find 
that r i/2=0.7" (corrected for seeing) and hence M^^IO 10 ' 3 
solar masses for Q0302— 223. 

In the case of Q1009— 0026, we know that the system is 
rotating and we can use the enclosed mass to determine the 
dynamical mass (Epinat et al. 2009): 



M d yn = Vmax ri/2 I G 



(2) 



where Vmax is the maximum velocity and rw2 is as be- 
fore. Estimating ri/2 from a curve of growth is more difficult 
in this case because of the proximity of the galaxy to the 
edge of SINFONTs field of view. As a consequence, we mea- 
sured Khmhm, the half width at halfmax, from a Gaussian 
fit to the ID profile extracted along the major kinematic 
axis. We find it to be Rh ro h m =0.64" , corresponding to rj/2 
of ~0.8" using the empirical relation ri/2 = l-25 x Khmhm for 
disks and under a seeing of 1" (Bouche et al. 2007b; Bouche 
in preparation). We therefore find Md yn =10 10 ' 9 solar masses 
for Q1009-0026 (see Table 2). 

Mass of Gas: In order to estimate the mass of gas in 
these objects, we start from the observed H-a surface bright- 
ness and compute gas surface brightness using an inverse 
"Schmidt-Kennicutt" relation (Bouche et al. 2007b, Finkel- 
stein et al. 2009): 



E gas [M e /pc 2 ] = 1.6 x 10" 



rg/s/kpc 2 ] 



(3) 



Gas masses inferred in this way are subject to signif- 
icant uncertainties beyond those associated with the mea- 
surements of fluxes and sizes and the 0.3 dex uncertainty 



SED Fit 




Wavelength (microns) 

Figure 2. SED fit of emission line galaxy to the broad-band 
magnitudes from HST/WFPC2 (Le Brun et al. 1997) and ground- 
based near-infrared observations (Chen & Lanzetta 2003) based 
on "Le Phare" photometric-redshift code (Ilbert et al. 2009). The 
spectroscopic redshift and star formation rate derived from our 
SINFONI spectra is used as an input to the code, thus allowing 
one to constrain the stellar mass of the object with relative high 
confidence, M»=10 9 ' 5 Mq and the age, about 2.5 Gyr. 



introduced by the scatter in the Schmidt-Kennicutt law it- 
self. It is not yet known whether the Schmidt-Kennicutt law 
in this form holds at high redshift. In addition, the calcula- 
tion is based on integrating the gas surface density over the 
visible region, but more gas could lie beyond this limit as in 
the case of the Milky Way. In the galaxy in the field of quasar 
Q0302— 223, we find a gas mass M 9as = 10 91 solar masses and 
we find 10 9 ' 2 in the other system. In comparison with the 
dynamical masses M^,, above, these results indicate a low 
gas fraction in the objects. In addition, both objects show 
clear exponential light profiles indicative of disks. 

Halo Mass: Thanks to our kinematical data based on 
the H-a emission widths (reflecting the velocity dispersion 
of the warm gas), we are able to estimate the mass of the halo 
in which the system towards Q1009— 0026 resides, assuming 
a spherical virialised collapse model (Mo & White 2002): 



M halo = O.Lff^G- 1 ^ -^! + z\ 



i.SyS 

1 ?"/) (7 X 



(4) 



using inclination-corrected Vmax values. We find 
M^ a i o =10 12,6 solar masses. This halo mass is compara- 
ble with the one from the Milky Way: 1.9±f 10 12 M 
(Wilkinson et al. 1999). 



5 SPECTRAL ENERGY DISTRIBUTION FIT 

Stellar Mass: In the case of the DLA towards Q0302-223, 
we use the broad-band magnitudes from HST/WFPC2 (Le 
Brun et al. 1997) and ground-based near-infrared obser- 
vations (Chen & Lanzetta 2003) which cover the Balmer 
break of the object to constrain the age of the stellar pop- 
ulation in the galaxy with a Spectral Energy Distribution 
(SED) fit to the integrated light of the galaxies. We used the 
photometric-redshift code "Le Phare" (Ilbert et al. 2009), 
with a single burst of star formation, templates from Bruzual 
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Table 2. Kinematic properties and mass estimates of the two N(H i) absorbers detected (see text for details of the computations). 
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0.88 


0.19 


0.7 


0.13 


10 10.3 


10 1.9 
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10 9.5 


Q1009-0026 


0.60 


1.45 


0.5 


0.31 


10 10.9 


10 2.2 


10 9.2 


10 12.6 





Note: The inclination is the main source of uncertainties and is estimated to be around 30%. 
a : The higher-resolution HST/WFPC2 data from Le Brun et al. (1997) clearly shows that the object is sub- 
divided into two sub-components, consistent with the elongated shape seen in the SINFONI data presented 
here. In this table, however, the object is treated as only one. 



& Chariot (2003) stellar population models, a Calzetti ex- 
tinction law and a Charbrier (2003) initial mass function. 
The spectroscopic redshift and star formation rate derived 
from our SINFONI spectra are used as an input to the code, 
thus allowing one to constrain the stellar mass of the object 
with relatively high confidence. We deduce M* = 10 9 ' 5 M 
for an object that would be about 2.5 Gyr old. An illustra- 
tion of the template which best fit the observed magnitudes 
is provided in Figure 2. 

This, in turn, allow us to put constraints on the baryonic 
mass fraction in this galaxy: 

fbaryons = (M gas + M t )/Mdyn (5) 

We derive a baryonic fraction of the order of 20%. Sim- 
ilarly, we can estimate the gas fraction is this object: 

H = M ga3 /(M gaa + M„) (6) 

We derive a gas fraction of l/3 rd . Such gas fractions 
are in the low range of the typical values derived in z~2-3 
galaxies by others (e.g. Erb et al. 2006, Law et al. 2009). 
When comparing these various mass estimates, we see that 
these systems have little room for molecular gas, which is 
consistent with the low star formation rates derived. 



6 COMPARISON WITH ABSORPTION 
PROPERTIES 

What makes our targets unique in comparison with other 
galaxies at z~l observed with Integral Field Unit (e.g. 
Forster-Schreiber et al. 2009, Law et al. 2009) is the wealth 
of information about the absorption properties of the gas 
available from high-resolution Echelle spectroscopy (see also 
Kacprzak et al. 2010). This allows us to compare the distri- 
bution of velocities of the neutral gas (absorption) and the 
ionised gas (emission). In the case of the absorber towards 
Q0302-223, Pettini et al. (2000) have fit the absorption 
profile with two main groups of components, separated by 
Av =36 km/s. They also indicate additional weaker compo- 
nents, at v=35 and 121 km/s relative to z a b s =1.009, which 
are visible in the stronger Fe II lines, n the case of the 
sub-DLA towards Q1009-0026, Meiring et al. (2006) have 
reported complicated velocity structure with seven compo- 
nents spanning Av ~334 km/s. This is considerably larger 
than the velocity dispersion measured at the position where 
it flattens out in the velocity map, a disk =60-70 km/s. This 



may support previous claims that the width of the absorp- 
tion systems cannot be solely explained by gravitational mo- 
tions but by the presence of outflows (Bouche et al. 2007a). 



7 CONCLUSION 

In conclusion, the observational set-up of SINFONI has 
demonstrated the power of integral field spectroscopy for de- 
riving a number of emission properties for quasar absorbers, 
a type of high-redshift galaxies that have been difficult to 
identify in the past. Detailed dynamical properties of these 
galaxies with known gas characteristics could be derived. 
These new tools are now available to systematically study 
these objects. We find that the two absorbers studied here 
have kinematical properties similar to other galaxies stud- 
ied in the same way at these and others redshifts (Epinat 
et al. 2009, Forster-Schreiber et al. 2006, Genzel et al. 2006, 
Forster-Schreiber et al. 2009). 

We find that galaxy associated with the DLA identified 
towards Q0302— 223 shows little signs of rotation and sig- 
nificant amounts of dispersions. Moreover, the light profile 
and morphology provide further evidence that this galaxy is 
dispersion dominated, albeit it already has a disky morphol- 
ogy. This is indicative of a young object which is confirmed 
by results from a SED fit or could be due to the blend- 
ing of the two components detected with HST interacting 
together. On the contrary, the galaxy associated with the 
sub-DLA towards Q1009— 0026 has a morphology and kine- 
matics consistent with that of a disk, with a normal disper- 
sion profile peaking at the center and flattening out in the 
outer-parts. This object shows clear signatures of rotation 
with systematic velocity gradients which is not typical of 
local disk galaxies, but the systematic gradient still favors 
a spiral galaxy. It remains to be understood why such high 
N(H i) column densities are found at large impact parame- 
ters (up to b=39kpc) from galaxies. 

Overall, we conclude that of the two absorbers stud- 
ied, the less metal-rich absorber toward Q0302— 223 arises 
in a gas-rich system with H-a luminosity consistent with 
that of a disk, while the more metal-rich absorber toward 
Q1009— 0026 arises in a lower gas mass system but with 
higher total mass showing clear signs of rotation. With the 
limitations of small number statistics, our findings are con- 
sistent with the suggestion (Khare et al. 2007; Kulkarni et 
al. 2007) that the metal-rich sub-DLAs may arise in more 
massive galaxies compared to the DLAs, which are usu- 
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ally metal-poor by comparison. Our work illustrates that 
detailed studies of quasar absorbers can offer entirely new 
insights into our knowledge of the interaction between stars 
and the interstellar gas in galaxies. 
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